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Abstract 

We present an analytical stability theory for the onset of the Faraday insta- 
bility, applying over a wide frequency range between shallow water gravity 
and deep water capillary waves. For sufficiently thin fluid layers the surface is 
predicted to occur in harmonic rather than subharmonic resonance with the 

forcing. An experimental confirmation of this result is given. 
PACS: 47.20.Ma, 47.20.Gv, 47.15.Cb 
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The observation of standing waves at the surface of a fluid layer subjected to a vertical si- 
nusoidal vibration dates back to Faraday []J . For sufficiently strong forcing the plane surface 
undergoes an instabilty giving rise to ordered patterns of different orientational order . 
For purely sinusoidal forcing spatially periodic patterns of lines, squares |5|||] , hexagons and 
triangles |7j have been found, even quasiperiodic structures of 8-fold orientational order have 
been observed ||. Faraday already recognized that the response of the surface appears with 
twice the period of the forcing. The first theoretical investigation of the linear stability 
H showed that the problem can be reduced to a set of Mathieu oscillators (parametrically 
driven pendula). However, the analysis takes advantage of the potential flow approximation 
which strictly applies to ideal (=inviscid) fluids only. Viscous effects are usually treated by 



introducing a heuristic damping term in the Mathieu equations |T(J being proportional to 
the kinematic viscosity v. Its strength is estimated by evaluating the viscous energy dissi- 
pation in the bulk of the fluid. This approximation ignores viscous boundary layers which 
appear close to the borders of the fluid. Within these boundary layers the potential flow 
assumption breaks down giving rise to vortical flow, associated with additional damping. 
Kumar and Tuckerman JTTJ] were the first who performed a linear stability analysis based 
on the viscid hydrodynamic field equations. Their fully numerical analysis uses Hill's in- 
finite determinant jl2| and allows to approximate the critical quantities with any desired 
accuracy. They pointed out that the stability problem for viscous fluids cannot rigorously 
be reduced to a damped Mathieu oscillator. Beyer and Friedrich demonstrated that a 
systematic treatment of the viscosity gives rise to a memory term in the Mathieu equation 
leading to an integro-differential equation. The integral gives rise to additional damping 
scaling - unlike bulk damping - with z/ 3 / 2 . Since Beyer and Friedrich used idealized free 
slip boundary conditions they did not catch damping in the bottom boundary layer, which 
becomes crucial - as we shall see - for small filling depths. Recently Kumar JTJJ presented 
an approximative linear analysis based on a truncation of his numerical method [|TTJ. For 
small damping he obtained an expression for the neutral stability curve, which requires a 
subsequent minimization with respect to the wave number to extract the threshold. Beside 
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of being rather implicit this approach is not systematic since one does not know to what 
order of viscosity the obtained result is valid Jl2|]. Since damping strongly increases as the 



wavelength A = 2n/k of the pattern compares to the filling depth h, Kumar |14j] points out 
the possibility of observing the harmonic surface resonance rather than the subharmonic 
one. The necessary parameter combination is beyond of his analytical approach, that's why 
he provides an example by numerical means. 

In the present paper we give a detailed account of the different damping mechanisms 
and develop a perturbative analytical treatment of the linear stability problem which does 
not suffer from the earlier shortcomings. We provide expressions for the critical onset accel- 
eration a c and the critical wave number k c applying also to the case of shallow water waves 
(A ~ h). The perturbation analysis is performed for the subharmonic (S) as well as for 
the harmonic (H) instability, allowing an easy prediction of the bicritical situation. Finally 
we present a Faraday experiment giving the first observation of surface waves in harmonic 
resonance with the forcing. 

We consider a layer of an incompressible fluid of density p, depth h with a free upper 
surface in contact with air. The layer is subjected to a sinusoidal vertical vibration corre- 
sponding to a modulated gravitational acceleration g(t) = g + acos2Qt in the comoving 
frame of reference. Control parameters are the strength of the modulation a and the forcing 
frequency / = 2Q/(2tc). The free surface is initially flat at the vertical coordinate z = 0. As 
the forcing amplitude exceeds a critical thershold a c the surface is located at the position 
z = ((x,y,t), where x,y are the horizontal coordinates and t is the time. For fluids of 



viscosity v and depth h the linear dispersion relation |14| for free (i.e. a = 0) surface waves 
of the form ( oc exp i(kx — out) is 

= A(k,u) = 

2 ^ r(r 4 + 2r 2 + 5) coth (rkh) — (1 + 6r 2 + r 4 ) tanh (kh) 2 
111 r coth (rkh) — coth (kh) 



4r(r 2 + 1) tanh (kh) / (cosh (kh) cosh (rkh)) 



r coth (rkh) — coth (kh) 



(1) 
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with co>o(A:) = tanh (kh) [g k + (a/p)k 3 ]/Q 2 , e = uk 2 /\Q\ and r = \Jl+ iu/e. The square 
root with the positive real part is assumed. To shorten notation we have nondimensionalized 
time t and frequencies u, Uq by the control parameter Q. In the ideal fluid limit v —>■ Eq. [T| 
reduces to the well known gravity-capillary dispersion relation u 2 = uj 2 (k). The periodic 
forcing is introduced by replacing go in Uq by g(t). This parametric driving couples different 
temporal Fourier modes and the following linear evolution equation for the Fourier transform 
((u) = j exp (—iut) ((t) dt of the surface elevation results 

A(k, W )C( W ) + akta ^ kh ) [ i{u _ 2) + i{u + 2 )] = . (2) 

It is useful to account for the different length scales of relevance in the Faraday experiment 
(see Fig. jT]). Intrinsic are the wavelength A = 2ir/k of the surface pattern, the thickness of 



the viscous boundary layer 5 = J2v/\Q\, and the capillary length Ja/(pg ), where a is the 
surface tension. Dictated by the geometry are the filling depth h and the lateral dimension 
of the vessel L, which we ignore here by assuming L — ► oo. 

In order to make analytical progress we confine ourselves to the limit of low viscosity e 
1, telling that the depth of the viscous boundary layer is small compared to the wavelength. 
Furthermore we assume that h is at least a few times larger than the thickness of the viscous 



boundary layer, h/^2u/\Q\ ~ 3, giving coth(rkh) ~ 1 and 1/ cosh (rkh) ~ 0. Note that 
there is no restriction upon the relation between h and A. The two simplifications made are 
not very restrictive and include almost all recent experiments. Even in the measurements 
of Edwards and Fauve with a high viscosity water-glycerol mixture e did not exceed a 
value of 0.4. We first expand Eq. [I] in powers of 1/r = 0(y/e) giving 

aii \ i i ■ /o i 2\ s x l 2 \fk + icu 3 

A(k, ijj) = uui - oo 2 + tuje(3 + coth 2 ) H h 

smh cosh 

e 3/ Ve + ^(-6tanh + coth + coth 3 ) + (3) 

where we have abbreviated coth (kh) by coth etc. This formulation is particularly useful for 
a physical interpretation. Transformation of Eq. 13] into real space yields a damped Mathieu 



oscillator with integral contributions |L5 
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= + e (3 + coth 2 )((t) + 

,1/2 



2 a k tanh 
+ — ^ — cos ( 2t ) 



C(t) 



£ 



7rsinh cosh 



/' G(t-r)(£ + 9 T ) 2 C(r)rfr 

J — oo 



-6 tanh + coth +coth 3 , /9 /"* 



e 3/2 / G(t-r)(£ + a T )C(r)rfr, 

J — oo 

(4) 

where G(t) = exp (— et)j\fl. Beside the usual (bulk) damping (oc Q the two integrals 
also contribute to the dissipation. This nonlocal temporal behavior is analogous to the 
propagation of thermal waves into a medium which is heated at its boundary. Here, the 
moving surface emitts velocity waves into the interior of the fluid giving rise to dissipation 
which depends on the history of C(^)- The first memory integral scales like Oiy 1 ! 2 ) and 
governs the expansion in the shallow water limit kh = 0(1). It is associated with damping in 
the bottom boundary layer and dies out exponentially (l/(sinh kh cosh kh)) as kh increases. 
The second integral (c.f. Ref. 0]) survives for h — > oo. It scales with u 3 ^ 2 and is related to 
the dissipation within the surface boundary layer. 

To proceed with the linear stability analysis we reexpand the square roots in Eq. |3| giving 

A(k,u) = -oo 2 +X(k,u) +Uq, (5) 

where all viscous contributions up to 0(e 3 ^ 2 ) are collected in 

X{k,uj) = &(X)+i$S(X) = 

1^13/2 M 1 / 2 

-e 1/2 1 1 h£ 3/ — ^(-15 tanh +5 coth +2coth 3 ) + 

V2sinhcosh 2y2 



1 1/2 

-\/2sinh cosh 2\/2 



£ 1/2 7' + e |u;|(3 + coth 2 ) + £ 3/2 ^^(-15 tanh + 5 coth + 2 coth 3 ) 



• (6) 



We first investigate the subharmonic resonance of the surface elevation by introducing 

CH = a x 5(u - 1) + a 2 5(uj + 1) + Ci + ••• (7) 

in Eq. ^|. Assuming a small detuning uo 2 — 1 = 0(X) the solvability condition for (j = O(X) 
yields the neutral stability curve 
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X(k, 1) + (wjj - 1)1 [*(*, -1) + K 2 - 1)1 = (^^J ■ (8) 



The minimum of the forcing amplitude a with respect to Uq ||T6[ relates the onset of the 
subharmonic response to the imaginary part of X 

of ~ — coth (fc/i) Of [X(k s , 1)] , (9) 
ks 

while the critical wave number k$ results from the dispersion relation corrected by the real 
part 

u 2 (k s )~l-M[X{k s ,l)) . (10) 
For a typical experimental situation Fig. compares our analytical formula Eq. with 



the full numerical computation according to the method of Kumar and Tuckerman []TT 
Excellent agreement is achieved over a wide frequency range between shallow water gravity 
and deep water capillary waves. The three contributions in 3(X) (Eq. |]) are respectively 
related to damping in the bottom boundary layer, the bulk, and the surface boundary 
layer, which become important at low, intermediate, and higher forcing frequencies. Even 
the sharp reincrease of the onset amplitude at low frequencies (shallow water limit kh ~ 
1) is well reproduced. Consistence of the perturbation expansion up to 0(e 3 ^ 2 ) requires 
l/(sinh cosh) < 0(e 1//2 ), which reflects a low frequency limit for the validity of Eqs. H,10. 



At high excitation frequencies Q the validity is limited by the size of the expansion coefficient 
e. 

A similar perturbation expansion can be computed for the first harmonic stability tongue. 
We obtain for the onset amplitude 

40 2 



a 



(H) 



k H 

with the critical wave number ku determined by 



coth(k H h)J$[X(k H ,2)} (11) 



u 2 (k H ) ~ 4 + jjs [X(k H , 2)] - 3? [X(k H , 2)] . (12) 
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A comparison with the exact numerical stability analysis (Fig. |2|) shows that the agreement is 
worse than for the subharmonic response. The sharp increase of at low frequencies leads 
to an intersection with a)* 1 ' and thus to a situation where the harmonic response preempts 
the subharmonic one. This phenomenon, however, is difficult to observe experimentally: 
The corresponding fi-window is rather narrow as it is also limited from below by the second 
subharmonic resonance (see Ref. [14|). An additional problem is a technical one and applies 



to most of the commercial shaker systems being in use for Faraday experiments: For the 
considered excitation frequencies (/ ~ 10Hz) it is usually the maximum peak elevation 
(but not the maximum force) which prevents the apparatus from reaching the threshold 
amplitude a^ H '. 

In order to observe the harmonic surface resonance we set up a Faraday experiment. We 
used a cylindrical aluminum container of radius R = 45mm and depth h = 5mm. Between 
R = 35mm and the outer edge the depth continuously decreases to zero. This provides 
increasing damping in order to suppress meniscus waves. Working fluid is a silicone oil 
(Dow Corning 200) with a viscosity of 8.9mm 2 / s, surface tension of 19.8* 10~ 3 iV/m, and 
density of 0.929g /cm 3 at a temperature T = 30° C. We use an electromagnetic shaker 
V400 (LDS) with a maximum force of 98iV and an indicated maximum peak elevation 
of 8 mm. Waveform generation as well as data aquisition is performed by a PC. The 
acceleration is measured by a piezoelectric device (Bruel and Kjaer 4393). For pattern 
visualization the container is lighted from the top by a concentric ring (20cm in diameter) of 
50 high intensity light emitting diods (LED) allowing for a stroboscopic illumination. A CCD 
camera located in the middle of the ring observes the pattern from the top. The LED's are 
synchronized to the forcing and triggered by either the same or twice the excitation period. 
Duty cycle as well as the relative phase of the illumination with respect to the forcing can be 
controlled externally. This technique allows a clear distiction between the subharmonic and 
harmonic surface resonance. As explained earlier bottom damping plays a crucial role for the 
appearance of the harmonic Faraday instability; thus a depth to wavelength ratio of kh ~ 1 
is required. For practical purposes filling depths much less than a millimeter are unsuitable 
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and difficult to control, we used h ~ 0.8mm. For this parameter combination the stability 
theory predicts the harmonic-subharmonic bicriticality, a[ s ^ = a{ H \ at / ~ 9Hz, which is 
close to the low frequency limit of our apparatus. Fig. |3]a presents a photograph of a pattern 
in harmonic resonance with the forcing. Both, the vessel and the surface were oscillating 
synchronously with 9Hz. To confirm our observation we show in Fig. |3|b a subharmonic 
surface pattern, oscillating again with 9Hz but excited with / = 18Hz. As expected, the 
wavelength in Fig. ^|a and b are the same. Finally Fig. |||c depicts a subharmonic pattern 
driven at / = 10Hz exhibiting a considerably larger wavelength. All patterns are achieved 
by slowly increasing the forcing amplitude a beyond the threshold while keeping / constant. 
The small aspect ratio of our container (diameter to wavelength ratio) did not allow the 
observation of nice ordered structures. Thus redoing the experiment with a larger (and 
thus havier) container would indeed be desirable, however, exceeds the specifications of our 
shaker [|17[ . 



In summary we have presented an analytical theory for the onset of the Faraday insta- 
bility. The analysis is based on the low viscosity approximation and assumes a filling depth 
larger than the thickness of the viscous boundary layer. Almost all recent experiments are 
covered by these assumptions. Particularly interesting is the case of shallow water waves 
(A ~ h) for which the harmonic instability preempts the subharmonic one. This theoretical 
prediction is confirmed by an experiment. 
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FIGURES 

FIG. 1. Length scales of relevance: Wavelength of the pattern A = 2ir/k, thickness of the 
viscous boundary layers 5 = y / 2z7/TT, depth of the layer h, and lateral extension of the container L. 
Energy dissipation in the various regions scales with different powers of the viscosity v (exponents 
are indicated in paranthesis). Waves with a wavelength larger (smaller) than the capillary length 
pgo/cr (not shown in the figure) are called gravity (capillary) waves. 

FIG. 2. The critical onset amplitude for the subharmonic (S) and harmonic (H) Faraday in- 
stability as a function of the forcing frequency / = Q/n. Solid lines correspond to the analytical 
results (Eqs.9 and 11), dashed lines are obtained by an exact numerical treatment. At low forcing 
frequencies the harmonic instability preempts the subharmonic one. Parameters: p = 0.934<7/cm 3 , 
cr = 20.1 * 10 _3 7V/m, v = 10mm 2 /s, h = 1mm. 

FIG. 3. Surface patterns as observed in the experiment, (a): Driving force and surface oscillate 
synchronously at / = 9Hz (harmonic response). ; (b) Subharmonic surface response with 9Hz at a 
forcing with / = 18Hz; (c) Subharmonic surface oscillation with 5Hz at a forcing with / = 10Hz. 




g.l: An analytical stability theory for Faraday waves 
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Fig. 2: An analytical stability theory for Faraday waves ... 

by Muller et al. 



1 1 1 1 r~7 — 1 I 1 1 1 1 i 1 1 1 1 i / 1 




_i i i i i i i i i_ 

50 150 250 



forcing frequency f = Q/% [Hz] 



Fig. 3: An analytical stability theory for Faraday waves 
by Muller et al. 



